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Abstract 



We establish sohton-Uke asymptotics for finite energy solutions to the Schrodinger 
equation coupled to a nonrelativistic classical particle. Any solution with initial 
state close to the solitary manifold, converges to a sum of traveling wave and out- 
going free wave. The convergence holds in global energy norm. The proof uses 
t:J- . spectral theory and the symplectic projection onto solitary manifold in the Hilbert 

I phase space. 
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1 Introduction 

We continue the study of coupled systems of wave fields and particles. In [10] the Klein- 
Gordon equation coupled to a relativistic particle has been considered. Here we extend 
the result to the Schrodinger equation coupled to a nonrelativistic particle. We prove the 
long time convergence to the sum of a soliton and dispersive wave. The convergence holds 
in global energy norm for finite energy solution with initial state close to the solitary 
manifold. 

We consider the Schrodinger wave function ip^x) in coupled to a nonrelativistic 
particle with position q and momentum p, governed by 



iip{x, t) = —Aip{x, t) + m^ijj{x, t) + p{x — q{t)) 
Qit) t)Vp{x - q{t)) + t)V^{x - q{t))\ dx, 

where m > 0. Denote tpi = Reip, ip2 = lnn-ip, pi = Rep, p2 = Imp. 
Then the system (1.1) becomes 

ipi{x,t) = -Aip2{x,t) + m^il)2{x,t) + P2{x - q{t)), 
^2{x,t) = Aipi{x,t) - m^ijji{x,t) - pi{x - q{t)), 

(iit) = / {Mx,t)Vpi{x - q{t)) + ij2{x,t)Vp2{x - q{t)))dx. 



xeR-" (1.1) 



xeR^ (1.2) 



This is a Hamilton system with the Hamilton functional 

7^(V'i,V2,g,g) = ^ J (|VV'i(x)|' + |VV'2(x)|' + m>i(x)|' + m>2(x)|') dx 

+ J (^Jpi{x)pi{x - q) +^2{x)p2{x - q)Jdx+^\q\'^. (1.3) 

We consider the Cauchy problem for the Hamilton system (1.2) which we write as 

Y{t)^F{Y{t)), teR; Y{0)^Yo. (1.4) 

HereF(t) = {i^i{t),ip2{t),q{t),p{t)), p{t) := q{t), Yq = (^oi, ^02, go,Po), and all derivatives 
are understood in the sense of distributions. Below we always deal with column vectors 
but often write them as row vectors. The system (1.2) is translation-invariant and admits 
soliton solutions 

YaA^) = i'^vi{x -vt - a), ijjv2ix - vt - a),vt + a,v), (1.5) 
for all a,v &R^ with \v\ < 2m. The states Sa,v '■= ^a,i.(0) form the solitary manifold 

S {Sa,v ■.a,ve K^ 1^1 < 2m}. (1.6) 
Our main result is the soliton asymptotics of type 

ijj{x,t) ijjyj^lx - v±t - a±) + Wo{t)ij^^, t — »■ ±00, (1.7) 

for solutions to (1.1) with initial data close to the solitary manifold S. Here ip^^ ~ 
'^v±i + ^V^D±2) WQ{t) is the dynamical group of the free Schrodinger equation, if^^ are the 
corresponding asymptotic scattering states, and the asymptotics hold in the global energy 
norm, i.e. in the norm of the Sobolev space H^{R^). For the particle trajectory we prove 
that 

q{t) ^ v±, q{t) v±t + a±, t — > ±00. (1.8) 



The results are established under the following conditions on the complex valued charge 
distributions p: 

{l + \x\fp, (l + |x|)^Vp, (l + |a:|)^VVpeL2(R=^), (1.9) 

with some /5 > 3/2. We require that all "modes" of the wave field are coupled to the 
particle, this is formalized by the Wiener condition 

p{k) = (27r)-^/2 j e''=^p(a;)cix ^ for all k . (1.10) 

It is an analogue of the Fermi Golden Rule: the coupling term p[x — q) is not orthogonal 
to the eigenfunctions e^^^ of the continuous spectrum of the linear part of the equation 
(cf. [4, 21, 22, 23]). 

Similar results were proved for the first time by Buslaev and Pcrelman [2, 3] for ID 
translation invariant Schrodinger equation, and extended by Cuccagna [6] for nD case, 
n > 3. In [10] the Klein-Gordon equation coupled to a particle, is considered. 

For the proofs of the asymptotics (1.7) and (1.8), we develop the approach [10] based 
on the Buslaev and Perelman methods [2, 3]: the symplectic orthogonal decomposition of 
the dynamics near the solitary manifold, the time decay for the linearized equation, etc. 
Our problem differs from [10] in the following aspects: 

i) Speed of propagation for the Schrodinger equation is infinite, and the solitons exist 
only for the velocities \v\ < 2m. 

ii) We consider nonspherically symmetric coupled function p{x). In this case we need 
additional arguments for the absence of embedded eigenvalues in the continuous spectrum. 

iii) We also consider the coupling function p{x) without compact support. Respec- 
tively, for the proof of the time decay for the linearized equation, we use the Jensen-Kato 
results [14, 15] and the Agmon weighted norms [1]. 

Remark 1.1. The term w? in Schrodinger equation appears automatically in the non- 
relativistic limit of the Klein- Gordon equation, and traditionally is removed by a gauge 
transformation. We keep the term to provide the existence of the nonzero solitons. 



2 Main Results 



2.1 Existence of Dynamics 

To formulate our results precisely, we need some definitions. We introduce a suitable phase 
space for the Cauchy problem corresponding to (1.2) and (1.3). Let = L^, and be 
the Sobolev space = {ip e L"^ : jV^'l G L^} with the norm = ||V^||l2 + ||^||l2. 

Let us introduce also the weighted Sobolev spaces H'^, s = 0, 1, o; e M with the norms 
||^||,,„:= ||(l + |x|)-^||^.. 

Definition 2.1. i) The phase space £ is the real Hilbert space (B (B R^ (B of 
states Y = {ipi, 1P2, q,p) with the finite norm 

||>'IU = ||^l||f/l + 11^21^1 + + 

a) £a is the space if^ © © © with the norm 

\\Y\\a = II YWs^ = ||V'l||l,a + ||V'2||l,a + |?| + N- 

iii) 8^ is space H"^ ® H"^ ® M? with the norm 

ll^llf = IIV'iIIh^ + ||V'2||i/2 + \q\ + bl- 



For ipj e we have 



^ ■llp.-lli^ < + ii^j, pA- - q)) < ^^ml^ + llMl^- (2.1) 



2m2"^^"^^ 2 '"^'"^ ^ 2 '"^'"^ 2' 

Therefore £^ is the space of finite energy states. The Hamihon functional Ti is continuous 
on the space S and the lower bound in (2.1) implies that the energy (1.3) is bounded from 
below. 

The system (1.2) reads as the Hamilton system 



/ 1 \ 

-10 

1 

\ -1 y 



r = (V^i,V'2,g,p)e£, (2.2) 



where VH is the Frechet derivative of the Hamilton functional (1.3). 
Proposition 2.1. Let (1.9) hold. Then 

(i) For every Yq ^ £ the Cauchy problem (1-4) has a unique solution Y{t) e C{M.,S). 
(a) For every t eM, the map U{t) : Yq Y{t) is continuous on £. 
(Hi) The energy is conserved, i.e. 

n{Y{t)) =n{Yo), teR. (2.3) 

Proof. Step i) Let us fix an arbitrary 6 > and prove (i)-(iii) for Iq £ ^ such that 
11^0 lU < b and \t\ < e — e{h) for some sufficiently small e{h) > 0. Let us rewrite the 
Cauchy problem (1.4) us 

y(i) = Fi(y(i)) + F2(y(i)), teR-. y(o) = yo, (2.4) 

where Fi : Y ^ {{~^ + ^^)^2) (^ ~ 'n2^)\l/i, 0, 0). The Fourier transform provides the 
existence and uniqueness of solution Yi{t) e C(]R, £^) to the hnear problem (2.4) with 
F2 — 0. Let Ui{t) : Yq i— > Yi{t) be the corresponding strongly continuous group of 
bounded linear operators on £. Then (2.4) for Y{t) e C(M, £) is equivalent to 

t 

Y(t)^U,{t)Yo + j dsU,{t-s)F2{Y{s)), (2.5) 


because F2{Y{-)) G C{M.,S) in this case. The latter follows from a local Lipschitz conti- 
nuity of the map F2 in 8: for each b > there exist a x = x(6) > such that for all 
Y,Z e£ with \\Z\\e < 6, 

\\F2{Y) - F2{Z)\\£ < h\\Y - Z\\£. 

Therefore, by the contraction mapping principle, equation (2.5) has a unique local solution 
F(-) G C{[—£,e],£) with £ > depending only on b. 

Step a) We use now energy conservation to ensure the existence of a global solution and 
its continuity. First consider Yq G := © Then Y{t) G since 
Ui{t)Yo, F2{Y{t)) G £+ by (1.9). The energy conservation law follows by (2.2) and the 
the chain rule for the Frechet derivatives: 

f^n{Y{t)) = {Dn{Y{t)),Y{t)) = {Dn{Y{t)), JDn{Y{t))) = O, teR 



since the operator J is skew-symmetric by (2.2), and DH{Y{t)) e ® ® ® for 
Y{t) e £:+. The inequahty (2.1) imphes 

1 1 1 

n > ^WWl^ + ^llV'lli. + - —Ml^- 

Hence, by energy conservation, for \t\ < e 

1 777^ 1 1 

2llvV'lli. + —\mh + ^\p\' - —JpfL^ < n{Y{t)) = mvo). 

This imphes a priori estimate 

\\j/j\\hi + \p\<B for \t\ < e, (2.6) 

with B depending only on the norm ||io||£ of the initial data and on ||p||i,2. An arbitrary 
initial data Yq E S can be approximated by initial data from Sc- The corresponding 
solution exists due to representation (2.5) by contraction mapping principle, and then 
(2.6) follows by the limit transition. 

Step Hi) Properties (i)-(iii) for arbitrary t £ R now follow from the same properties for 
small \t\ and from a priori bound (2.6). □ 

2.2 Solitary Manifold and Main Result 

Let us compute the solitons (1.5). The substitution to (1.1) gives the following stationary 
equations 

-iv ■ Vipy{y) = (-A + m'^)ijy{y) + p{y) 

p = v, ^ = -J {Wviy)piy) + ^My)piy))dy ^^'^^ 

Then the first equation implies 

^My) ■= [-^ + m'' + iv v]My) = -p{y), y e R'. (2.8) 

For \v\ < 2m the operator A is an isomorphism i?''(R^) —>■ H^{M.^). Hence (1.9) implies 
that 

i^M = -A-'piv) e H^W"). (2.9) 

If V is given and 1^1 < 2m, then py can be found from the second equation of (2.7). 

The function t/jy can be computed by the Fourier transform. The sohton is given by 
the formula 

M^) = / \ \ 2.10 

47r J \x- y\ 

Further, in Appendix A, we prove that the last equation of of (2.7) holds. Hence, the 
soliton solution (1.5) exists and defined uniquely for any couple {a.v) with \v\ < 2m. Let 
us denote by := e R^ : 1^1 < 2m}, ipyi = Reipy, and ipy2 = Iraipy. 

Definition 2.2. A soliton state is S{a) :— {ipvi{x — b),ipv2{x — b),b,v), where a :— {b, v) 
with 6 e R^ and v eV. 

Obviously, the soliton solution admits the representation S{a{t)), where 

a{t) = {b{t),v{t)) = {vt + a,v). (2.11) 
Definition 2.3. A solitary manifold is the set S :— {S{a) : cr e E := R^ x y}. 



The main result of our paper is the following theorem. 



Theorem 2.1. Let (1.9), and the Wiener condition (1.10) hold. Let (3 > 3/2 be the 
number from (1.9), and Y{t) be the solution to the Cauchy problem (1.4) with the initial 
state Yq which is sufficiently close to the solitary manifold: 

po < 2m, do := dist£^{Yo, S) < 1. (2.12) 

Then the asymptotics hold for t — > ±oo, 

m =v± + Oi\t\-'), qit) = v^t + a± + Oi\t\-^/'y, (2.13) 



^{x, t) = ^^±{x - v±t - a±) + W^o(^)V'± + r±{x, t) (2.14) 

with 

\\r±mm^O{\t\-'/'). (2.15) 

It suffices to prove the asymptotics (2.13), (2.14) for t — > +oo since the system (1.2) 
is time reversible. 



3 Symplectic Projection onto Solitary Manifold 

Let us identify the tangent space to S, at every point, with £. Consider the symplectic 
form fl defined on £ by O J dipi{x) A dip2{x) dx + dq A dp, i.e. 

n{Y,,Y2) = {Y,,JY2), Yi,Y2eS, (3.1) 

where 

(Yi, Y2) := (-011, -012) + (-021, -022) + qiq2+PiP2 
and {'ipu,'4'i2) = j '4^ii{x)ipi2{x)dx etc. It is clear that the form Q is non-degenerate, i.e. 

n{Yi,Y2) = for every Y2 e £ =^ Yi = 0. 

Definition 3.1. i) Yi \ Y2 means that Yi e Y2 £ ^, o,nd Yi is symplectic orthogonal to 
Y2, I.e. fi(Yi,Y2) = 0. 

a) A projection operator P : S ^ £ is called symplectic orthogonal if Yi { Y2 for 
Yi e KerP andY2 G ImP. 

Let us consider the tangent space Ts[a)S to the manifold <S at a point S{a). The 
vectors Tj :— d(j.S{a), where d^. :— df,. and ^o-^+g := d^. with j — 1, 2, 3, form a basis in 
XjS. In detail. 



= ^ji'") ■= db^S{a) = ( -dji)yi{y) , -djijjy2{y) , ej , 
'Tj+3 ^ 'Tj+si'^) — dy.S{a)^{ d^.i/j^iiy) , 9^^.V'^2(?/) , 0, Cj 



J = 1,2, 3, (3.2) 



where y := x — b is the "moving frame coordinate", ei = (1,0,0) etc. Let us stress that 
the functions tj will be considered always as the functions of y, not of x. 
The formulas (2.10) and the conditions (1.9) imply that 

Tj{v)eSa, veV, j = l,...,6, Va</?. (3.3) 



Lemma 3.1. The matrix with the elements D,{Ti{v),Tj{v)) is non- degenerate for any v e 
V. 

Proof The elements are computed in Appendix B. As the result, the matrix ^l{Ti,Tj) 
has the form 

niv) (Q(r,,r,)),,=,,...,6 = ( f ) , (3.4) 

where the 3 x 3-matrix fl'^{v) equals 

n+(v)^K + E. (3.5) 
Here is a symmetric 3 x 3-matrix with the elements 

^ J (F + m2)2 - (A;^;)2 ^^"^^ 

where the "hat" denotes the Fourier transform (cf. (1.10)). The matrix K is the integral 
of the symmetric nonnegative definite matrix k <S> k = {kikj) with a nonnegative weight. 
(The last statement is true since k"^ + > \{kv)\ for \v\ < 2m, and j'^^i + i'0t,2p = 

\'4^v2\'^ — ii'ipvii^v2 ~ '4^v2'ipvi) — 0-) Hcncc, the matrix K is also nonnegative definite. 
Since the unite matrix E is positive definite, the matrix Q'^(v) is symmetric and positive 
definite, hence non-degenerate. Then the matrix Q{ri,rj) also is non-degenerate. ■ 
Let us introduce the translations : {ipi{-),ip2{-),(l,p) ^ (V'i(" ~ct))V'2(- ~ o,), Q + o,,p), 
a e R^. Note that the manifold S is invariant with respect to the translations. 

Definition 3.2. i) For any a G M andp < 2m denote by Sa{p) — {Y — {ipi,ip2, q,p) € 

■■ \p\ < P}- We set £{p) := So{p). 
a) For any v < 2m denote by S(l;) = {a = (6, f ) : 6 e R^, \v\ < v}. 

The next Lemma provide that in a small neighborhood of the soliton manifold S a 
"symplectic orthogonal projection'' onto S is well-defined. The proof is similar to the 
proof of the Lemma 3.4 in [10]. 

Lemma 3.2. Let (1.9) hold, a G M. Then 

i) there exists a neighborhood Oa{S) of S in and a map 11 : Oa{S) S such that 11 
is uniformly continuous in the metric of Sa on Oa{S) fl Saip) with p < 2m, 

UY = Y for Kg 5, and Y-S\rsS, where S = YIY. (3.7) 

a) Oa{S) is invariant with respect to the translations Ta, and 

UTaY = TaUY, for Y G Oa{S) and a G R^ 

Hi) For any p < 2m there exists a v < 2m s.t.IlY — S{a) with a G S(v) for Y G 

o^{S) n s^ip) 

iv) For any v < 2m there exists an ra{v) > s.t. S{a) + Z E Oa{S) if a E E(t;) and 

\\Z\\a < ra{v). 

We will call 11 a symplectic orthogonal projection onto S. 

CoroUeiry 3.1. The condition (2.12) implies that Yq — S -\- Zq where S — S{ao) — HYq, 
and 

\\Zoh < 1. (3.8) 



Proof Lemma 3.2 implies that nio = "S" is well defined for small do > 0. Furthermore, 

the condition (2.12) means that there exists a point Si ^ S such that \\Yq — SiWp = do- 
Hence, Yq, Si G Oj3{S) r\Sf3{p) with ap< 2m which does not depend on do for sufficiently 
small do- On the other hand, II^'i = 5*1, hence the uniform continuity of the map 11 
implies that H^i — S\\f) ^ as do 0. Therefore, finally, \\Zo\\p = \\Yo — S\\p < 
\\Yo - Si\\p + \\Si - S\\p <do + o(l) < 1 for small do- ■ 



4 Linearization on the Solitary Manifold 

Let us consider a solution to the system (1.2), and split it as the sum 



Y{t)^S{a{t)) + Z{t), 



(4.1) 



where a{t) = {b{t),v{t)) e E is an arbitrary smooth function of i e R. In detail, denote 
Y = {tpi, q,p) and Z = (*i, *2, Q, P)- Then (4.1) means that 



V^i(a;,t) = V'^(t)i(x-6(t)) + ^i(x-6(t),t), q{t) = b{t) + Q{t) 
M^,t) = A{t)2ix-b{t)) + ^2{x-b{t),t), pit) = v{t) + P{t) 



(4.2) 



Let us substitute (4.2) to (1.2), and linearize the equations in Z. Later we will choose 
S{a{f)) = IIY(t), i.e. Z(t) is symplectic orthogonal to Ts[a{t))S. 

Setting y = X — b{t) which is the "moving frame coordinate", we obtain from (4.2) 
and (1.2) that 



^1 = v-V^^l:^i{y)-b-V^pyi{y) + <^i{y,t)-b-V^i{y,t) 

= - AV'^s (y) + m'^iPM - ^My, t) + ^'*2(y, t) + P2{y- Q) 

ij2 = V- Vv^v2{y) - b ■ V^pv2{y) + ^2{y, t)-b- V*2(z/, t) 

= AV'^i(y) - m''iP^i{y) + A^i{y, t) - m''^i{y, t) - pi{y - Q) 

q = b+Q^v+P 

p = ^, + p = _(v(V;,,•(|/) + ^',(l/,^)),p,•(l/-g)). 



(4.3) 



Let us to extract linear terms in Q. First note that pj{y — Q) = Pj{y) — Q ■ ^ Pjiy) + 
Nj{Q), J = 1,2, where Nj{Q) = pj{y - Q) - Pj{y) + Q ■ Vpj{y). The condition (1.9) 
implies that for Nj{Q) the bound holds. 



\Nj{Q)\\o,p<C0{Q)Q', i = l,2. 



(4.4) 



uniformly in < Q for any fixed Q, where /3 is the parameter from Theorem 2.1. Using 
the equations (2.7), we obtain from (4.3) the following equations for the components of 



the vector Z{t): 

^i{y,t) = -A^2{y,t)+m^^2{y,t)+b-V^,{y,t)-Q-Vp2{y) 

+ (b-v)-VMy)-v-V^'ijj^i(y) + N2 

^2{y,t) = A^i{y,t)-m''^i{y,t) + b-V^2{y,t) + Q-Vpi{y) 

+ {b-v)-yijy2{y)-v-yvi^v2{y)-Ni 

Q{t) = p + {v-h) 

P{t) = {^lfj{y,t),Vpj{y)) + {V^,j{y),Q-Vpj{y))-v + N^{v,Z) 

where N^{v,Z) = -{ViP^j, Nj{Q)) + (V*,-,Q • Vpj) - {V^j,Nj{Q)). Clearly, N4v,Z) 
satisfies the following estimate 



(4.5) 



N4{v,Z)\ < Cp{p,v,Q) + ||^i||i,-/3|Q| + ||^2||i,-/3|g| 



(4.6) 



uniformly in \v\ < v and \Q\ < Q for any fixed v < 2m. We can write the equations (4.5) 

as 



Z{t) = A{t)Z{t) +T{t) + N{t), teR. 

Here the operator A{t) = A^,«,(t) depends on two parameters, v = v{t), and w 
and can be written in the form 











Q 




\ P ) 





( w-V -(A-m^) -Vp2- 

A - w • V Vpr 

^ 

V (^vpi) (•,vp2) {v^vj,-vpj) y 



\ / ^1 \ 

^2 

Q 



(4.7) 

-m 



(4.8) 



V p J 

Furthermore, T{t) = T^,u,{t) and N{t) = N{t,a,Z) in (4.7) stand for 



T 



( {W -V) ■ ViPyl - V ■ Vyipyl \ 

{W -v) ■ - V ■ Vvi^vl 

V — W 
—V 



N{a,Z) 



( N2{Z) \ 
-Nr{Z) 


V N^{v,Z) ) 



(4.9) 



where v — v{t), w — w{t), a — a{t) — {b{t),v{t)), and Z — Z{t). The estimates (4.4) and 
(4.6) imply that 



\\Nia,Z)y<Civ, 



(4.10) 



uniformly in cr G S(v) and < r_^(v) for any fixed v < 2m. 

Remark 4.1. i) The term A{t)Z{t) in the right hand side of the equation (4.7) is linear 
in Z{t)^ and N(t) is a high order term in Z{t). On the other hand, T{t) is a zero order 
term which does not vanish at Z{t) = since S{a{t)) generally is not a soliton solution if 
(2.11) does not hold (though S{a{t)) belongs to the solitary manifold), 
n) Formulas (3.2) and (4.9) imply: 



T{t)^-Y^[{w-v)iri + viri+3] 



(4.11) 



and hence T{t) G Ts(a(t))<S, t G K. This fact suggests an unstable character of the 
nonlinear dynamics along the solitary manifold. 



5 The Linearized Equation 

Here we collect some Hamiltonian and spectral properties of the generator (4.8) of the 
linearized equation. First, let us consider the linear equation 

X{t) = A^^^X{t), teR, V eV, w eM.^. (5.1) 

Lemma 5.1. (cf. [10]) i) For any v & V and w e the equation (5.1) can be written 
as the Hamilton system (cf. (2.2)), 

X{t) = JDn,,^{X{t)), t e R, (5.2) 

where D7iv,w is the Frechet derivative of the Hamilton functional 

n^,n,{X) = ^ J [|V*i|' + m2|*i|2 + |V*2r + m2|*2r](i2/ + J *2«^ • V*idy 

+ y"p.(?/)Q-V^,dy + ip2-^(g-V^,,(y),g-Vp,(|/)), X = i^,,^2,Q,P)e£, 

(5.3) 

a) Energy conservation law holds for the solutions X{t) e C^(M, £"*"), 

HvA^it)) = const, t e R. (5.4) 
Hi) The skew-symmetry relation holds, 

^{Av,w^l, ^2) = —^{^l,Av,w^2), ^1,^2 G ^- (5.5) 

iv) The operator Ay ^y, acts on the tangent vectors Tj{v) to the solitary manifold as follows, 

Ay^^[Tj{v)] ^ (w-v)- VT,-('t;), Ay^y,[Tj+3{v)] = (w - -y) • "s/Tj+siv) + Tj{v), j = 1, 2, 3. 

(5.6) 

We will apply Lemma 5.1 mainly to the operator corresponding to w — v. In 
that case the linearized equation has the following additional essential features. 

Lemma 5.2. Let us assume that w = v & V . Then 

i) The tangent vectors Tj{v) with j = 1,2,3 are eigenvectors, and Tjj^^iv) are root vectors 
of the operator Ay corresponding to zero eigenvalue, i.e. 

Ay,y[Tj{v)]^Q, Ay,y[Tj+3{v)\^Tj{v), i = l,2,3. (5.7) 

a) The Hamilton function (5.3) is nonnegative definite since 

'>ivAX) = \ j |A'/'(*i + ^^2)-A"'/'g-V(pi+^p2)rcia; + ^p2>0. (5.8) 

Here A is the operator (2.8) which is symmetric and nonnegative definite in L^(M^) for 
\v\ < 2m, and A^/^ is the nonnegative definite square root defined in the Fourier repre- 
sentation. 



Proof. The first statement follows from (5.6) with w — v. In order to prove ii) we rewrite 
the integral in (5.8) as follows: 

^{A'^^^i + ^^2) - A-'/'Q • V(pi + ^p2), A^/'(^i + ^^2) - A'^/^g . y ^ ^^^^^ 

= l(A(*i + i*2), *i + ?*2) - i^j, Q ■ Vpj) + ^(A-'g • V(pi + ip2), Q ■ V(pi + ip2)) 

(5.9) 

since the operator A^/^ is symmetric in L^(M^). Now all the terms of the expression (5.9) 
can be identified with the corresponding terms in (5.3) since 

^(A(*i + W2), *i + i*2) = + m'^ + iv- V](*i + ^^2), (*i + ^^2)) 

= ^([-A + m2]^i, ^1) + ^([-A + m^]^2, ^'2) + (^2, v ■ V^i) 

and A-\pi + ip2) = -{ipvi + i'4^v2) by (2.8) and (2.9). □ 

Remark 5.1. For a soliton solution of the system(1.2) we have 6 = v, v = 0, and hence 
T(t) = 0. Thus, the equation(5.1) is the linearization of the system (1.2) on a soliton 
solution. In fact, wc do not linearize (1.2) on a soliton solution, but on a trajectory 
S{a{t)) with a{i) being nonlinear in t. We will show later that T{t) is quadratic in Z{t) 
if we choose S{a{t)) to be the symplectic orthogonal projection of Y{t). Then (5.1) is 
again the linearization of (1.2). 



6 Symplectic Decomposition of the Dynamics 

Here we decompose the dynamics in two components: along the manifold S and in 
transversal directions. The equation (4.7) is obtained without any assumption on a{t) in 
(4.1). We are going to choose S{a{t)) :— ny(i), but then we need to know that 

Y{t)e0^p{S), teM, (6.1) 

It is true for t = by our main assumption (2.12) with sufficiently small do > 0. Then 
S{a{0)) = UY{0) and Z(0) = F(0) - S{a{0)) are well defined. We will prove below that 
(6.1) holds if do is sufficiently small. Let us choose an arbitrary such that |t'(0)| <v< 2m 
and let S = v — 1^(0)1. Denote by r_^(u) the positive numbers from Lemma 3.2 iv) which 
corresponds to a = —f3. Then S{(7) + Z G (9_/3(iS) if a = {b,v) with \v\ < v and 
< r^isiy). Note that ||Z(0)||_/3 < r^piy) if do is sufficiently small. Therefore, 
S{a{t)) = UY{t) and Z{t) = Y{t) - S{a{t)) are well defined for i > so small that 
1^1 < V and ||Z(i)||_^ < r^p{v). This is formahzed by the following standard definition. 

Definition 6.1. U is the "exit time", 

t^ = sup{t>0:\\Z{s)\\^f3<r_p{v), \v{s) - v{0)\ < S, 0<s<t}. (6.2) 

One of our main goals is to prove that = 00 if do is sufficiently small. This would 
follow if we show that 

\\Z{t)\\_p<r_i3{v)/2, \v{s)-v{Q)\<d/2, 0<t<t,. (6.3) 

Note that _ 

\Q{t)\<Q ■■=r-/3{v), 0<t<t.. (6.4) 

Now N{t) in (4.7) satisfies, by (4.10), the following estimate, 

\\N{t)\\^<C^{v)\\Z{t)tf„ 0<t<t.. (6.5) 



6.1 Longitudinal Dynamics: Modulation Equations 

Prom now on we fix the decomposition Y{t) — S{a{t)) + Z{t) for < i < by setting 
S{a{t)) — TlY{t) wfiicfi is equivalent to tfie symplectic ortfiogonality condition of type 
(3.7), 

Z{t)\Ts^,^t))S, 0<t<U. (6.6) 

This aUows us to simphfy drastically the asymptotic analysis of the dynamical equations 
(4.7) for the transversal component Z{t). As the first step, we derive the longitudinal 
dynamics, i.e. the "modulation equations" for the parameters cr(t). Let us derive a 
system of ordinary differential equations for the vector cr(t). For this purpose, let us write 
(6.6) in the form 

n(Z(t),r,(t)) = 0, J = 1,...,6, 0<t<U, (6.7) 

where the vectors Tj{t) — Tj{a{t)) span the tangent space Ts(^a{t))^- Note that a{t) — 
{b{t),v{t)), where 

\v{t)\ <v <2m, 0<t<t„ (6.8) 

by Lemma 3.2 iii). It would be convenient for us to use some other parameters {c,v) 

instead of cr = {b,v), where c{t) — b{t) — / v{T)dT and 

Jo 

c{t) = b{t) -v{t) =w{t) -v{t), 0<t<U. (6.9) 

We do not need an explicit form of the equations for (c, v) but the following statement, 
which can be proved similar to the Lemma 6.2 in [10]. 

Lemma 6.1. LetY{t) be a solution to the Cauchy problem (1-4), o.nd (4-1), (6-7) hold. 
Then {c{t),v{t)) satisfies the equation 

where 

N{a,Z)^0{\\Z\\%) (6.11) 

uniformly m cr e E(v). 

6.2 Decay for the Transversal Dynamics 

In Section 11 we will show that our main Theorem 2.1 can be derived from the following 
time decay of the transversal component Z{t): 

Proposition 6.1. Let all conditions of Theorem 2.1 hold. Then t^ = oo, and 

mt)U<^^,, t>0. (6.12) 

We will derive (6.12) in Sections 7-11 from our equation (4.7) for the transversal 
component Z{t). This equation can be specified using Lemma 6.1. Indeed, the lemma 
implies that 

\\Tit)\\p<Civ)\\Zit)f_,, 0<t<t,, (6.13) 

by (4.9) since w — v = c. Thus (4.7) becomes the equation 

Z{t) ^ A{t)Z{t) + N{t), 0<t<t^, (6.14) 



where A{t) — and N{t) :— T{t) + N{t) satisfies the estimate 

\\N{t)y<C\\Z{t)f_^, Q<t<U. (6.15) 

In all remaining part of our paper we will analyze mainly the basic equation (6.14) to 
establish the decay (6.12). Wc are going to derive the decay using the bound (6.15) and 

the orthogonahty condition (6.6). 

First, we reduce the problem to the analysis of the frozen linear equation, 

X(t)^AiX(t), teR, (6.16) 

where Ai is the operator defined by (4.8) with vi = v{ti) and a fixed ti G [0,it*). 

Then we can apply wellknown methods of scattering theory and then estimate the error 
by the method of majorants. 

Note, that even for the frozen equation (6.16), the decay of type (6.12) for all solutions 
does not hold without the orthogonality condition of type (6.6). Namely, by (5.7) the 
equation (6.16) admits the secular solutions 

3 3 

X{t) = ^C,T,{v) + J2DArAv)t + r,M] (6-17) 
1 1 

which arise also by differentiation of the soliton (1.5) in the parameters a and v in the 
moving coordinate y — x — Vit. Hence, we have to take into account the orthogonality 
condition (6.6) in order to avoid the secular solutions. For this purpose we will apply 
the corresponding symplectic orthogonal projection which kills the "runaway solutions" 
(6.17). 

Remark 6.1. The solution (6.17) lies in the tangent space Ts^ai)^ with ai — (&i, vi) (for 
an arbitrary bi G R) that suggests an unstable character of the nonlinear dynamics along 
the solitary manifold (cf. Remark 4.1 iii)). 

Definition 6.2. i) For v E V, denote by n„ the symplectic orthogonal projection of £ 

onto the tangent space Ts(a)S, and P^, = I — 11^,. 

ii) Denote by = P^S the space symplectic orthogonal to Ts{a)S with a = {b,v) (for an 
arbitrary b eM.). 

Note that by the hnearity, 

n„z = Yl ^A^M^Mniv), z), ze£, (6.18) 

with some smooth coefficients Tlji{v). Hence, the projector 11^,,, in the variable y = x — b, 
docs not depend on b, and this explains the choice of the subindex in 11^, and P^,. 
Now we have the symplectic orthogonal decomposition 

£^rs^,)S + Z,, a^{b,v), (6.19) 

and the symplectic orthogonality (6.6) can be written in the following equivalent forms, 

n^(t)Z(t)=0, P^^t)Z{t)^Z{t), 0<t<t,. (6.20) 

Remark 6.2. The tangent space Ts{a)S is invariant under the operator A^^y by Lemma 
5.2 i), hence the space is also invariant by (5.5): Ay^^Z e Z^ for sufficiently smooth 
Z eZy. 



In Sections 12-18 we will prove the following proposition which will be one of the main 

ingredients for proving (6.12). Let us consider the Cauchy problem for the equation (6.16) 
with A — Ay^y for a fixed v & V. Recall that the parameter /9 > 3/2 is also fixed. 

Proposition 6.2. Let the conditions (1.9) and (1.10) hold, \v\ <v< 2m, and Xq e £. 

Then 

i) The equation (6.16), with A = A^^^, admits the unique solution c'^^Xq := X{t) G 

C{M.,S) with the initial condition X{0) = Xq. 

a) For Xq & Zy n Sp, the solution X{t) has the following decay. 



\\^^'^o\\-p<j^^^^\\Xoh, tew. (6.21) 
7 Frozen Transversal Dynamics 

Now let us fix an arbitrary ti G [0,t^,), and rewrite the equation (6.14) in a "frozen form" 

Z{t) = AiZ{t) + {A{t) - Ai)Z{t) + N{t), 0<t<t„ (7.1) 
where Ai = Ay^ti),v{ti) 

/ [w{t) - v{ti)] • V \ 

[w{t) - v{t,)] ■ V 

A{t)-A,- Q Q 0- 

\ (V(^.(t),--V'.(ti).),Vp,) 0/ 

The next trick is important since it allows us to kill the "bad terms" [w{t)—v{ti)] ■ V in 
the operator A{t) — Ai. 

Definition 7.1. Let us change the variables {y,t) i— > {yi,t) — {y + di{t),t), where 



di{t) := [ {w{s) -v{ti))ds, 0<t<ti. 



(7.2) 



Next define 

Z,{t) ■.= {^,{y,-d^it),t),^2{yl-d^it),t),Qit),Pit)). (7.3) 
Then we obtain the final form of the "frozen equation" for the transversal dynamics 

Zi{t)^AiZi{t)+Bi{t)Zi{t)+Ni{t), 0<t<ti, (7.4) 
where Ni{t) — N{t) expressed in terms of y = yi — di{t), and 

/ \ 





V (V(^,(t),--z^„(i,),-),Vp,-) J 

Lemma 7.1. (see [10]) For (^i, ^2, Q, -P) G with any a < (3 the following estimate 
holds: 

||(*i(2/i - til), *2(yi - cii), g, P)|U < *2, Q, P)Ui + Mi|)i"i , e (7.5) 

CorollEiry 7.1. The following hounds hold for < t < ti 

\\Ni{t)h < \\z{t)r_p{i + \d,{t)\r , \\B^{t)z,{t)y < c\\z{t)uf^^ \\z{T)\\%dT . 

(7.6) 



8 Integral Inequality 



The equation (7.4) can be written in the integral form: 

= e^i*Zi(0) + [ e^'^^-''>[BiZi{s) + Ni{s)]ds, Q<t<ti 
Jo 



(8.1) 



Now we apply the symplectic orthogonal projection Pi := Pv(ti) to both sides of (8.1). 
The space Zi :— Pi£ is invariant with respect to e^^* by Proposition 6.2 ii) (cf. also 
Remark 6.2). Therefore Pi commutes with the group e^^* and applying (6.21) we obtain 
that 

IIP 7 mil ^r ^lPiZMh , r t \\Pi[BiZ,{s) + m{s)]\\0 ds 

||PlZi(t)||_,<C ^^^^^3/^ +Cj^ (1 + It -.1)3/2 ■ 

The operator Pi = I — Hi is continuous in £^ by (6.18). Hence, using 7.6 we obtain that 



+ C(di(t)) 



0(l + |t-s|)3/2 



\Z(s) 



\ziT)r_pdT+\\zis)r_p 



ds, 0<t< ti. 

(8.2) 



where di{t) := supo<s<t Mi(s)|. 
Definition 8.1. t'^ is the exit time 

< = sup{t e [0,t,) : rfi(s) < 1, < s < t}. 
Now (8.2) imphes that for ti < t'^ 
C 



.3) 



|Pi^i(^)ll-^<^j:p^ll^(o)||, 



+ Ci 



r L 

io (1 + 1^- 



; 1)3/2 



\Zis) 



\Z{T)\\%dT+\\Z{s)\\% 



ds, 0<t< ti. (8.4) 



9 Symplectic Orthogonality 



Finally, wc arc going to change PiZi(t) by Z{t) in the left hand side of (8.4). We will 
prove that it is possible using again that do <^ 1 in (2.12). For the justification we reduce 
further the exit time. First, we introduce the "majorant" 



m{t) sup {l + sf/^\\Z{s)\\^p , t e [Q,t,). 

se[o,t] 

Let us denote by £ a fixed positive number which we will specify below. 
Definition 9.1. t" is the exit time 

t'i = sup{t e [0, <) : m{s) <e, 0<s<t}. 



(9.1) 



(9.2) 



The following important bound (9.3) allows us to change the norm of PiZi(t) in the 
left hand side of (8.4) by the norm of Z{t). 



Lemma 9.1. (cf.flOj) For sufficiently small e > 0, we have for ti < t" 

\\Z{t)\\.p < C||PiZi(t)||_;j, 0<t<h, (9.3) 

where C depends only on p and v. 

Proof. Since < 1 for t <ti < t" < t'^ then by Lemma 7.1 it suffices to prove that 

||Zi(i)||_;3<2||PiZi(i)||_^, Q<t<h. (9.4) 
Recall that PiZi(t) = Zi{t) — n^(t^)Zi(t). Then estimate (9.4) will follow from 

\\U,^t,)Zi{t)\\-P < l\\Zim-P: 0<t<h. (9.5) 
Symplectic orthogonahty (6.20) implies 

n,(t),iZi(t) = 0, te[o,ti], (9.6) 

where IIy(^t),iZi(t) is n^(t)Z(t) expressed in terms of the variable yi = y + di{t). Hence, 
(9.5) follows from (9.6) if the difference Hv{ti) — ^v{t),i is small uniformly in t, i.e. 

||n,(t,)-n,(t),i|| < 1/2, o<t<h. (9.7) 

It remains to justify (9.7) for small enough £ > 0. Formula (6.18) implies the following 
relation 

n,(i),iZi(t) = J2^jMt)h,Mt)Mri,Mt)),z,{t)), (9.8) 

where Tj,i(f(t)) are the vectors Tj{v{t)) expressed in the variables yi. Since |(ii(t)| < 1 
and Vtj are smooth and fast decaying at infinity functions. Lemma 7.1 implies that 

tj,Mt)) - rj(v(t))y < C\di(t)f, 0<t<h (9.9) 

for all j = 1, 2, . . . , 6. Furthermore, 

^jHt)) - ^jHti)) = v{s) ■ V^Tj{v{s))ds, 

and therefore 

\\Tj{v{t)) -Tj{v{ti))\\f3 <C \v{s)\ds, 0<t<ti. (9.10) 

Similarly, 

nji{v{t))-nji{v{h))\ = \J^' v{s)-v^nji{v{s))ds\ <c J^' \v{s)\ds, o<t<ti, (9.11) 

since \VvIlji{v{s))\ is uniformly bounded by (6.8). Hence, the bounds (9.7) will follow 
from (6.18), (9.8) and (9.9)-(9.11) if we establish that \di{t)\ and the integral in the right 
hand side of (9.10) can be made as small as we please by choosing £ > small enough. 
To estimate di{t), we note that 

w{s) - v{ti) = w{s) - v{s) + v{s) - v{ti) = c{s) + v{T)dT (9.12) 



by (6.9). Hence, (7.2), Lemma 6.1 and the definition (9.1) imply that 
\di{t)\ = I {w{s) - viti))ds\ < ' (^\cis)\ + ^ ' \vir)\dr^ ds 

< Cm\u) ^ ' (y ^ilsf + / ' (^^""^^3 ) ds < Cm\u) < Ce\ 0<t<h (9.13) 
since ti < t". Similarly, 

/ \v{s)\ds<Cm^{ti) ^<Ce^, 0<t<ti. (9.14) 

□ 



10 Decay of Transversal Component 

Here we prove Proposition 6.1. 

Step i) We fix e > and t'l = t'l{e) for which Lemma 9.1 holds. Then the bound of type 
(8.4) holds with ||PiZi(t) ||_^ in the left hand side replaced by ||Z(t)||_^ : 

H 1 r /■*! 1 

ds, 0<t<ti (10.1) 



for ti < t'^. This implies an integral inequality for the majorant ■m{t) defined in (9.1). 
Namely, multiplying both sides of (10.1) by (1 + t)^/^, and taking the supremum in 
t e [0, ti], we get 



m{ti) < C||Z(0)||^+C sup 



(1 + tf/^ 



m{s) rn?{T)dT w?{s) 



+ 



{l + sf/^J, (l + r)3 (l + s)3 



ds 



teloMJo {l + \t-s\r/^ 
for ^1 < t". Taking into account that m(t) is a monotone increasing function, we get 

m{h)<C\\Z{0)y + C[m^{h)+m\h)]I{ti), h<tl (10.2) 

where 

(i + i)3/2 r 1 /■*! 



r {i + tf^ r 1 dr i 

^ tIZJo {i+\t-s\r/n{i+sr/'Js (i+r)3 + (i+.)3^ 



ds < I < 00. 



te[o,ti] . 

Therefore, (10.2) becomes 

m{h) < C\\Z{0)y + C7[m^{h)+m''{h)], h < C (10.3) 
This inequality implies that m{ti) is bounded for ti < t", and moreover, 

m{h)<C,\\Z{0)\\p, h<t:, (10.4) 
since m(0) = ||Z(0)||^ is sufficiently small by (3.8). 

Step ii) The constant Ci in the estimate (10.4) does not depend on i*, t'^ and t" by Lemma 



9.1. We choose do in (2.12) so small that ||^(0)||^ < e/{2Ci). It is possible due to (3.8). 
Then the estimate (10.4) implies that t'l = t'^ and therefore (10.4) holds for all ti < t'^. 
Then the bound (9.13) holds for all t < t'^. We choose e so small that the right hand side 
in (9.13) does not exceed one. Then t[ = t*. Therefore, (10.4) holds for all ti < t*, hence 
the first inequality in (6.3) also holds if ||Z(0)||/3 is sufficiently small by (9.1) and (9.14). 
Finally, this imphes that — oo, hence also — t'^ — oo and (10.4) holds for all ti > 
if do is small enough. ■ 



11 Soliton Asymptotics 

Here we prove our main Theorem 2.1 under the assumption that the decay (6.12) holds. 
First we will prove the asymptotics (1.8) for the vector components, and afterwards the 
asymptotics (1.1) for the fields. 

Asymptotics for the vector components From (4.3) we have q — b + Q, and from 
(6.14), (6.15), (4.8) it follows that Q^P + 0(||Z||2^). Thus, 

q = b + Q = v{t) + c{t) + P{t) + C(||Z||%). (11.1) 
The equation (6.10) and the estimates (6.11), (6.12) imply that 

m\ + m\<^f!^, t>o. (11.2) 

Therefore, c{t) = c+ + 0{t-'^) and v{t) = v+ + ^(r^), t oo. Since |P| < the 
estimate (6.12), and (11.2), (11.1) imply that 

q{t)=v+ + 0{t-'/-'). (11.3) 

Similarly, 

b{t)^c{t)+ / v{s)ds^v+t + a+ + 0{t-^), (11.4) 
Jo 

hence the second part of (1.8) follows: 

q(t) = b(t) + Q{t) = v+t + a+ + O(r'), (11.5) 

since Q{t) = 0{t-^/^) by (6.12). 

Asymptotics for the fields We apply the approach developed in [12], see also [10]. For 
the field part of the solution, iplx, t) = ipi{x, t) + iip2{x, t) let us define the accompanying 
soliton field as ipv{t){x — q{t)), where we define now y{t) — q{t), cf. (11.1). Then for the 
difference z{x, t) — ip{x, t) — '0v(t)(2^ — ?(^)) we obtain the equation 

iz{x,t) = (-A + m^)z{x,t) - iv ■ Vyi>y{t){x - q{t)). 

Then ^ 

z{t) = Wo{t)z{Q) - [ Wo{t - s)[w(s) • VvV'v(.)(- - q{s))]ds. (11.6) 
Jo 

To obtain the asymptotics (2.14) it suffices to prove that z{t) = Wo{t)ijj_^_ + r+(t) with 
some e and ||r+(t)||j^i = 0{t^^^^). This is equivalent to 



Wo(-t)z(t)^^^ + r'4t), 



(11.7) 



where ||r^(t)||ifi = 0{t~^^'^) since W^{t) is a unitary group in the Sobolev space by the 
energy conservation for the free Schrodinger equation. Finally, (11-7) holds since (11-6) 
implies that 

Wo{-t)z{t) = z{0)- [ Wo{-s)f{s)ds, /(s) = ^v(s)•Vv^v(.)(•-g(s)), 

where the integral in the right hand side of converges in the Hilbert space with the rate 
0{t-^/^). The latter holds since \\Wo{-s)f{s)\\Hi = Ois-^/"^) by the unitarity of Wo(-s) 

and the decay rate ||/(s)||j^i = 0(s~^/^). Let us prove this rate of decay. It suffices to 
prove that |v(s)| = 0{s~^^^), or equivalently |p(s)| = Substitute (4.2) to the 

last equation of (1.2) and obtain 

Pit) = J [4^vit)jix - b{t)) + ^Hjix - b{t),t)] V pj{x-h{t)-Q{t))dx = J i^v(t)jiy)^Pj{y)dy 

+ J M)jiy) [Vpj{y-Q{t))-Vpj{y)]dy + J ^j{y,t)Vpj{y-Q{t))dy. 

The first integral in the right hand side is zero by the stationary equations (2.7). The 
second integral is 0{t~^^'^), since Q{t) = (9(t~^/^), and by the conditions (1.9) on p. 
Finally, the third integral is 0{t~^^'^) by the estimate (6.12). The proof is complete. ■ 



12 Decay for the Linearized Dynamics 

In remaining sections we prove Proposition 6.2 in order to complete the proof of the main 
result (Theorem 2.1). Here we discuss our general strategy of the proof of the Proposition. 
We apply the Fourier-Laplace transform 

POD 

X{X)= e-^'X{t)dt, ReA>0 (12.1) 
Jo 

to (6.16). According to Proposition 6.2, we expect that the solution X(t) is bounded in 
the norm || • ||_^. Then the integral (12.1) converges and is analytic for ReA > 0. We 
will write A and v instead of Ai and Vi in all remaining part of the paper. After the 
Fourier-Laplace transform (6.16) becomes 

XX(X) ^AX(X)+Xo, ReA>0. (12.2) 

Let us stress that (12.2) is equivalent to the Cauchy problem for the functions X{t) e 
C6([0, oo); £_^). Hence the solution X{t) is given by 

X{X) ^ -{A- Xy^Xo, ReA>0 (12.3) 

if the resolvent R{X) — {A — X)~^ exists for Re A > 0. 

Let us comment on our following strategy in proving the decay (6.12). First, we will 
construct the resolvent -R(A) for Re A > and prove that it is a continuous operator in 
Then X{X) G S^p and is an analytic function for ReA > 0. Second, we have to 
justify that there exist a (unique) function X{t) G C([0, cxd); £^_^) satisfying (12.1). 

The analyticity of -^(A) and Paley- Wiener arguments (see [16]) should provide the 
existence of a - valued distribution X{t), t eW, with a support in [0, oo). Formally, 

X{t) = 7^ /" e^^X{iu; + Q)du;, t e M. (12.4) 



However, to check the continuity of X{t) for t > 0, we need additionally a bound for 
X(iu} + 0) at large \uj\. Finally, for the time decay of X(t), we need an additional infor- 
mation on the smoothness and decay of X{iu! + 0). More precisely, we should prove that 
the function X{iu! + 0) 

i) is smooth outside a; = and uj = ±fj,, where /x = fj.{v) > 0, 

ii) decays in a certain sense as — > oo. 

iii) admits the Puiscux expansion at a; = ±/x. 

iv) is analytic at = if Xq G := P^S and Xq G Sf^. 

Then the decay (6.12) would follow from the Fourier-Laplace representation (12.4). 

We will check the properties of type i)-iv) only for the last two components Q{X) and 
P{X) of the vector X{X) = (^i(A), §2(A), Q{X), P{X)). The properties provide the decay 
(6.12) for the vector components Q{t) and P{t) of the solution X(t). Then for the field 
components ^i(a;, i) and ^2(2^, ^) we will use wellknown properties of free Schrodinger 
equation. 



13 Constructing the Resolvent 



Here we construct the resolvent as a bounded operator in for Re A > 0. We will write 
(^i(l/),^2(2/),Q,-P) instead of (^i(|/. A), ^^2(2/, A), (^(A), P(A)) to simplify the notations. 
Then (12.2) reads 



(^-A) 



^2 
Q 

\ P J 



^02 

V ^0 y 



It is the system of equations 

V ■ V<^!i{y) - (A - m')-^2{y) - Q ■ Vp2 - X^i{y) = -^foiiv) 
(A - m2)*i(y) + V ■ V*2(y) + Q • Vpi - A^'2(y) = -*02(y) 

P-XQ = -Qo 

-{V^j{y),pj(y)) + {V'ilj,j(y),Q- Vpj(y)) - AP = -Pq 
Step i) Let us study the first two equations. In Fourier space they become 
-{ikv + X)i>i{k) + (/c^ + m^)i>2{k) = -^oi{k) - zQkp2, 
-{e + m^)ifi{k) - {ikv + X)if2{k) = -*02(A:) + iQkp^. 
Let us invert the matrix of the system and obtain 



y e 



;i3.i) 



(13.2) 



f —{ikv + X) k'^ + w? \ ^ r/., ,x2 /,2 2\2-]-i( — {ikv + X) —{k'^ + m'^)\ 
\-{k'' + w?) -{ikv + X) J - il^/^^+Aj +^ ) J A;2 + ^2 _^ikv + X) )■ 

Taking the inverse Fourier transform we obtain the corresponding fundamental solution 



where 

(13.4) 

Note that denominator in RHS (13.4) does not vanish for ReA > 0, /c G M.^. Moreover it 
does not vanish for Re A > 0, A; e for sufficiently small | Im/cj. Therefore gx{y) decays 
exponentially by the Paley- Wiener arguments. Let us compute the entries of matrix Gx 
explicitly: 

G\\y) = Gf{y)=F-'-- — (13.5) 



1-1 



l/2i l/2i 



''^'^Kk'^ + m'^ — kv + iX k? + -\- kv — iX) 8i7r|y| 8i7r|y| 

_ / 1/2 1/2 \ (.-'<+\y\-i^y ^-^-\y\+qy 

-^h — i-i/ (to. o t , • a TO. o . T • A 1 I I 



'^^^ V A;2 _|- 177,2 _ _l_ ^_)^ /^2 _|_ ^2 _|_ _ 87r|y| 87r|y 



where 



This implies 



^2 



x± = -y/m^ - — ± iA, Rex± > 0. (13.6) 



Lemma 13.1. i) The operator G\ with the integral kernel Gx{y — y'), is continuous 
operator //^(M^) © H^[R^) ^ H'^{m^) ® H'^{R^) for Re A > 0. 

a) The formulas (13.5) and (13.6) imply that for every fixed y , the matrix function Gx{y), 
Re A > 0, admits an analytic continuation in X to the Riemann surface of the algebraic 
function ^/Jjfi^rX? with the branching points X — ±i/x, where // := — ^. 

Thus, from (13.2) and (13.3) we obtain the convolution representation 

= _G'ii*^Q^_Gi2^^^2-(Gf *Vpi)-g+(Gii*Vp2)-g, (13.7) 
*2 = G"^ * ^01 - Gl' * ^02 - (Gl' * Vpi) • Q - (Gi^ * Vp2) • Q- 

Step a) Now we proceed to the last two equations (13.1): 

-AQ + P = -Qo, {ViPvj, Q ■ Vpj) - (V^',-, pj) - AP = -Po. (13.8) 

Let us rewrite equations (13.7) as = j{Q) + ^j(^oi) ^02)5 where 

*i(*oi, *02) = -G"^ * *oi - G^^ * *02, *i(Q) = (-Gi^ * Vpi + * Vp2) • Q, 

*2(*oi, ^02) = * *oi - G]^ * *02, *2((5) = -{G]^ * Vpi + * Vp2) • Q. 

Then (V^j,Pj) = (V^'j(Q), p^) + (V^'j(^'oi, ^02), Pi), and the last equation (13.8) be- 
comes 

(VV'.i, Q ■ Vp,) - (Vvl/,(g), p,) - AP = -Po + (V^,(^oi, ^02), Pi) =: -Po- 
First we compute the term 

(VV'^,j, Q ■ Vpj) = "^{S/iJvj, QApj) = '^{'^i'vj, diPj)Qi. 
ij ij 



Applying the Fourier transform Fy_,k, we have by the Parseval identity and (19.3) that 
^{dii^vj^dipj) = ^(-ikiij^j^-ikipj) (13.9) 



— {k+m)pi+ikvp2 , —ikvpi — {k+m)p 



-/ 



/i;j/i;/(|(/i;2 + m2)(|pi|2+ \p2\^) + i{kv){piP2 - p2pi)^dk 
{k^ + m^y-{kvy 



-: —La. 



As the result, (Vipvj,Q ■ Vpj) = —LQ, where L is the 3x3 matrix with the matrix 
elements Lu. Now let us compute the term — (V^j((5), Pj) = {^j{Q), Vpj). One has 

(*,(g), diPj) = Yl {(-^'' * ^IP^ + * dlP2: diP,) - {G'^ * dip, + G'^ * diP2, diP2))Ql 

I 

I 

and again by the Parseval identity we have 

Hu{\): = {-G^^*dipi + G]^*dip2,d,pi)-{G]^*dipi + G^^*dip2,diP2) (13.10) 
= ([(/c^ + m^)pi - {ikv + \)p2]g\ki, hpi) + {[{ikv + A)pi + {k^ + m^)p2]g\ku hh) 
kiki({k'^ + m^)(|pip + IpaH + {ikv + A)(pi^2 - P2pi)^dk 



I 



The matrix H is well defined for Re A > since the denominator does not vanish. As the 
result, — (V^j(Q),pj) = i/Q, where H is the matrix with matrix elements Hn. Finally 
the equations (13.8) become 

-^l^) ( ? ) = ( ?I ) ■ ^(^) = ( L -H(X) aI ) • (".11) 

Assume for a moment that the matrix A4{X) is invertible (later we will prove this). Then 
we obtain 

(p)=A^-'(A)(^?), ReA>0. (13.12) 

Finally, formula (13.12) and formulas (13.7), where Q is expressed from (13.12), give the 
expression of the resolvent -R(A) = (^4 — A)~^, Re A > 



Lemma 13.2. The matrix function A1(A) (respectively, Ai ^(A)j; ReA > admits an 
analytic (respectively meromorphic) continuation to the Riemann surface of the function 
^Jp? + A2, A e C. 

Proof. The analj^ic continuation of M.{\), exists by Lemma 13.1 ii) and the convolution 
expressions in (13.10) by (1.9). The inverse matrix is then meromorphic since it exists for 
large ReA. The latter follows from (13.11) since H{X) ^ 0, ReA ^ oo, by (13.10). □ 



14 Analyticity in the Half-Plane 



Here we prove the following 



Proposition 14.1. The operator-valued function R{X) : S ^ S is analytic for ReA > 0. 

Proof. It is sufficient to prove that the operator A — X : £ ^ £ has a bounded inverse 
operator for ReA > 0. Let us recall, that A = Ay^y where \v\ < 2m. 

Step i) Let us prove that Ker (A — A) = for ReA > 0. Indeed, let us assume that 
Xx = (*ai,*a2,Qa,-Pa) e ^ satisfies (A - X)Xx = 0, that is Xx is a solution to (13.1) 
with \E'oi = ^02 = and Qo = Po = 0. We have to prove that Xx = 0. 

First let us check that Pa = 0. Indeed, the trajectory X := Xxe^^ G C{M.,£) is the 
solution to the equation X = AX that is (5.1) with w = v. Then T-ly^y{X{t)) grows 
exponentially by (5.8). This growth contradicts to the conservation of Tiy^y, which follows 
from Lemma 5.1 ii) since X{t) e C^(R, £^+). The latter inclusion follows from Lemma 
13.1 since (^ai, ^A2) satisfies the equations (13.7) with ^oi = *02 = and Q — Qx. 

Now AQa = -Pa = by the third equation of (13.1), hence Qa = since A 7^ 0. Finally, 
^^Ai = 0, ^'a2 = by the equations (13.7) with Q = Qx = 0. 

Step ii) Let us represent A — X — Aq -\-T, where 



/ vW-X -(A - m^) 



Ao-- 



\ 








V- 





A 





-A 




\ 




-a; 



/ 














-•Vp2 

•Vpi 





\ 


E 
J 



The operator T is finite-dimensional, and the operator A^^ is bounded in £ by Lemma 
13.1. Finally, A — X = Aq{I + A^^T), where A^^T is a compact operator. Since we know 
that Ker (7 + A^^T) = 0, the operator (7 + A^^T) is invertible by Fredholm theory. □ 



Corollary 14.1. The matrix M{X) of (13.11) is invertible for Re A > 0. 



15 Regularity on the Imaginary Axis 

First, let us describe the continuous spectrum of the operator A — A^^^ on the imaginary 
axis. By definition, the continuous spectrum corresponds to cu G M, such that the resolvent 
i?(ic(j + 0) is not a bounded operator in £. By the formulas (13.7), this is the case when the 
Green function Gx{y — y') loses the exponential decay. Thus, iur belongs to the continuous 
spectrum if 

\uj\ > fj, — rr? — v'^ /A. 

By Lemma 13.2, the limit matrix 

exists, and its entries are continuous functions of a; G R, smooth for < // and \uj\> ^. 
Recall that the point A = belongs to the discrete spectrum of the operator A by Lemma 
5.2 i), hence M.{iuj + 0) (probably) also is not invertible at a; = 0. 

Proposition 15.1. Let p satisfy the condition{1.9) and the Wiener condition (1.10), and 
\v\ < 2m. Then the limit matrix A4{iu! + 0) is invertible for ou ^ 0, a; G R. 



Proof We will consider separately three cases < \uj\ < fx, u — fx, and > fj,. We 
can assume that v = {\v\,0,0). Let us denote F{lj) := 



o = + |P2p, b — i{piP2 — P2Pi). Then the entries of the matrix F become 

1 1 



F ■ 



Ma 



+ b 





v\ki+ cu 


M2 - 


-( 


V 





/Iv'} k-idk ( 



M2- (|t;|A;i+u;)2 - {\v\kxY 

\v\kx 



+ 



L + i7(ia; + 0), M = m^ + k'^, 

(15.2) 



1 



1 



1 



+ b 



1 



M-\v\ki-uj M+\v\ki+u M-\v\ki M + \v\ki 
1 1 1 



M-\v\ki-u M+\v\ki + uj M-\v\ki M + \v\ki 
Since a is even, and b is odd we obtain that 

+00 

dk2dk^ J kikj dki 



afi + bf2 



(15.3) 



where 



/i: = 



/2: = 



1 



1 



1 



+ 



1 



M-\v\ki-u; M+\v\ki+u M + \v\ki - u M-\v\ki+u; 
2 2 



(15.4) 



M-\v\ki M+\v\ki 
1 1 



+ 



1 



M-\v\ki-u; M+\v\ki + u; M-\v\ki+u; M + \v\ki - u 

9 9 

+ 



M-|t;|A;i M+|t;|A;i' 
Then by (15.1) 



det M{iu!) = det 



/ iu -1 \ 

iuj 0-10 

io; -1 

-Fn -Fi2 -Fi3 icj 

-Fi2 -F22 -F23 iij 

\ -Fi3 -F23 -F33 y 

3 



= -a,''-io' 



UJ 



i<j 



Fn 


F12 


Fi3 


F12 


F22 


F23 


Fi3 


F23 


F33 



(15.5) 



since Fij = Fji. 

I. First, let us consider the case < \u\ < fx. Then the invertibihty of A4{iu;) follows 
from 

Lemma 15.1. For < < fx, the matrix F is positive definite. 

Proof. First, let us note that all denominators in (15.4) are positive for |ci;| < fx — rri^ — 
v'^/A, \v\ < 2m. Indeed, 

.,2 .,2 

2 I _2 , , „. /n\2 , _2 



(m' + ey - (u; + \v\k^y = {{k - vj^y + m 



u){{k^vl2y ^m^ - — + a;) > 



Second, fi>f2>0if\v\ < 2m and < |a;| < /x. This is proved in Appendix C. 
Finally, the Wiener condition implies 

a±b=\pi{k)Tip2ik)\'^ > 0, Vfc G Ml (15.6) 

Therefore afi + 6/2 > and (15.3) is the integral of the symmetric nonnegative definite 
matrix k®k — {kikj) with a positive weight. Hence, the matrix F is positive definite. □ 

II. uj = ±/x. Let us consider for example uj = n = — In this case (13.10) reads: 

kikj{Ma - {kv + iJi)h)dk 



{{k, - M)2 + + kl) ((A;i + M)2 + ki + kl + 2/.) ■ 



Now the integrand has a unique singular point. The singularity is integrable, hence 
detM.{iu}) also is negative by the representations (15.5). Hence, the matrix is also 

invertible. 

III. \uj\> 11. Here we apply an other arguments. Now the invertibility of M.{iuj) follows 
from (15.5) by the following lemma (cf. [10]) 

Lemma 15.2. // (1.10) holds and uj > /i (uj < —n), then the matrix ImF(a;) is negative 
(positive) definite. 

Proof. We consider the case u > n (the case u < can be treated similarly). Let us 
calculate the imaginary part of F^j. Since F^j = Hij{iuj + 0) — L^j and L^j is real, we will 
consider only Hij{iu + 0). For £ > we have 

, . f kikj{Ma + {kv + u — ie)b)dk 1 f kikj{a + h)dk . 

H,^{iuj + e) = j M2_^kv + uj- lef =27 M-kv-uj^ie ^^^'^^ 
If kk,{a-h)dk ^ . . 
2jM + kv + u-t£ ' ' 

It suffices to consider only the first summand in (15.7), since the second summand is real 
for £ = 0. Consider the denominator 

Deik) ^ k^ + rn^ - kv - u; + ie. 

Note that Do{k) = on an ellipsoid T^, where 

Then the Plemelj formula for C^-functions implies that 

ImHUzuj + 0) = - J / (15.8) 

where dS is the element of the surface area. Hence, the matrix Im H^{iuj + 0) is negative 
definite by (15.6). □ 



Now let us prove that the limit matrix M.{iuj + 0) is invertible. Recall that 

M{iuj + Q) = 



luE -E 
-F(ia; + 0) iuE 



Then the equation 



M{iu; + 0) ( p ) = 



becomes 

iujQ-P = 0, -FQ + iujP = 0. (15.9) 

Then {F + uj'^)Q = 0, which imphes Q = and then P = since the matrix ImF is 
negative definite for uo > ix. This completes the proofs of the Proposition 15.1. 

Corollary 15.1. Proposition 15.1 implies that the matrix Ai~^{iuj) is smooth m a; e M 
outside three points a; = 0, ±/x. 



16 Singular Spectral Points 

Let us recall that the formula (13.12) expresses the Fourier-Laplace transforms Q{\), -P(A). 
Hence, the components are given by the Fourier integral 

{%)=ije-'M-\^^0)[<f,y (16.1) 

if it converges in the sense of distributions. The Corollary 15.1 alone is not sufficient for 
the proof of the convergence and decay of the integral. Namely, we need an additional 
information about a regularity of the matrix M~^{iuj) at its singular points a; = 0, ±//, 
and some bounds at — > oo. We will analyze all the points separately. 

I. First we consider the points ±/i. 

Lemma 16.1. The matrix A4^^{iuj) admits the following Puiseux expansion in a neigh- 
borhood of±ii: there exists an e±> s.t. 

oo 

M'\iu;) ^J2^k(^^l^)''^^^ \oo^l^\<e±, u eR. (16.2) 

fe=0 

Proof. It suffices to prove a similar expansion for Ai{iu!). Then (16.2) holds also for 
M~^{iu;), since the matrices M{±ifJ,) are invertible. The asymptotics for M.{iu;) holds 
by the convolution representation (13.10) 

H,j{X) = -{Gf * dip, + Gl' * dip2,d,p,) - {Gl' * dip, + G'^ * dip2,d,p2). (16.3) 

since G^^ admit the corresponding Puiseux expansions by the formula (13.5). □ 

II. Second, we study the asymptotic behavior of A^~^(A) at infinity. Let us recall that 
M.~^{X) was originally defined for Re A > 0, and admits a meromorphic continuation to 

the Riemann surface of the function ^ m"^ — ^ + iX (see Lemma 13.2). 

Lemma 16.2. There exist a matrix Rq and a matrix-function Ri{w), such that 

R 

M'^iiu) = — + Ri(u;), \u;\>ii+l, u eR, 

UJ 

where, for every A; = 0, 1, 2, .. 

|a^i?i(a;)| < \uj\>p+l, ueR. (16.4) 



Proof. The structure (15.1) of the matrix M.{iuj) provides that it suffices to prove the 
following estimate for the elements of the matrix H{iuj) :— H{iu; + 0): 



\dtHjj{tu;)\ < 



\UJ\ 



\uj\>l^+l, j = 1,2,3. 



;i6.5) 



Note, that 



G\'*f^ l(Dr^(A)/ - D-\X)f), G'^*f^ l{D^\X)f + D^\X)f), 



where 



Di{X) = -A + m^ - i-y • V + iA, D2{X) = -A + + • V - iA, Re A > 0, 



and Dj (A), j — 1,2 are bounded operators in 
follows from a more general bound 

\\^tDJ\^u + Q)f\\,.<^ 



The estimate (16.5) immediately 



\UJ\ 



|a;| > /X + 1 



;i6.6) 



which holds for a > 3/2. Namely, (16.5) follows by (1.9) from (16.6) applied to the 
functions f{y) = dipj{y) e Ll. 

The bound (16.6) is proved in [1, the bound (A.2')] (see also [15, Thm 8.1]). □ 

III. Finally, we consider the point cu — which is most singular. The point is an 
isolated pole of a finite degree by Lemma 13.2, hence the Laurent expansion holds. 



M~'^{iuj) ^'^MkUj-''-^ + n{uj), \uj\<eo. 



(16.7) 



k=0 



where Mk are 6x6 complex matrices, Eq > 0, and Ti-ico) is an analytic matrix- valued 
function for complex cu with < Eq. 



17 Time Decay of the Vector Components 

Here we prove the decay (6.12) for the components Q{t) and P{t). 

Lemma 17.1. (cf. [10]) Let Xq e Zy^p. Then Q{t), P{t) are continuous and 



m)\ + \p{t)\< 



C{p,v,do) 

(1 + 1^1)3/2' 



t > 0. 



(17.1) 



Proof. The expansions (16.2), (16.4) and (16.7) imply the convergence of the Fourier 
integral (16.1) in the sense of distributions to a continuous function of t > 0. Let us prove 
(17.1). First let us note that the condition Xq G implies that the whole trajectory 
X{t) lies in Zy^p. This follows from the invariance of the space Z^^/j under the generator 
Ay^y (cf. Remark 6.2). Note that for Xq not belonging to Zy^p the components Q{t) 
and P{t) may contain non-decaying terms which correspond to the singular point a; = 0. 
Indeed, we know that the linearized dynamics admits the secular solutions without decay, 
see (6.17). The formulas (3.2) give the corresponding components Qs{t) and Ps{t) of the 
secular solutions. 



Qs(t) 
Psit) 



EC.- 







;i7.2) 



We will show that the symplectic orthogonality condition leads to (17.1). Let us split the 
Fourier integral (16.1) into three terms using the partition of unity Ci((^) +C2(<^) +C3('^) — 
1, a; e M: 



= hit) + hit) + hit), (17.3) 
where the functions Cki^) G C°°(R) are supported by 



supp Ci 


c 


{a; e 1 


: £0/2 < \uj\ < 11 + 2} 


supp C2 


c 




:\uj\> 11+1} 


supp Cs 


c 


{a; e M 


■ kl < £0} 



(17.4) 



Then 

i) The function hit) decays, hke (1 + |i|)~^/^, by the Puiseux expansion (16.2). 

ii) The function hit) decays faster than any power of t due to Proposition 16.2. 

iii) Finally, the function hit) generally does not decay if n > in the Laurent expansion 
(16.7). Namely, the contribution of the analytic function TCiuj) decays faster than any 
power of t. On the other hand, the contribution of the Lorent series, 

( ^ '^^'^^^^^ |j ^^^^ " ( pI ) ^ ^ ^' ^^^-^^ 

is a polynomial function of i e R, of a degree < n, modulo functions decaying faster than 
any power of t. Let us note that the formula (17.2) gives an example of the polynomial 
functions appeared from (17.5). 

We have to show that the symplectic orthogonality condition eliminates the polynomial 
functions. Our main difficulty is that we do not know anything about the order n of the 
pole and the Lorent coefficients of the matrix A4~^(ia;) at a; = 0. 

Our crucial observation is the following: 

a) The components (17.2), of the secular solutions, form a linear space Cs of the dimen- 
sion dim £5 = 6. 

b) The polynomial functions from (17.5) belong to a linear space of the dimension 
dim Cl < 6, since (Qo, ^o) e K^- 

c) Cs C Cl since all the functions (17.2) admits the representation (17.5). The latter fol- 
lows from the fact that the secular solutions (6.17) can be reproduced by our calculations 
with the Laplace transform. 

Therefore, we conclude that 

Cl = Cs. (17.6) 

It remains to note that the secular solutions are forbidden since Xq e 2^v,/3- Hence, the 
polynomial terms in (17.5) vanish that implies the decay (17.1). 

More precisely, we know thatX(t) = P^Xit) for all t G M. On the other hand, the 
identity (17.6) implies that X(t) can be corrected by a secular solution Xsit) s.t. the 
corresponding components (5a(^) and -Pa(^), of the difference A(t) := Xit) — Xsit), 
decay. Hence, the components Qit) and Pit), of Xit) — P^X(t) = P^[X(t) — Xsit)], also 
decay. □ 



18 Time Decay of Fields 

Here we prove the decay of the field components ^^1(2;, t), ^^2(2^, t) corresponding to (6.12). 
The first two equations of (6.16) may be written as one equation: 

i^(t) = (-A + m^ + iv- V)* - Q{t) ■ Vp, (18.1) 

where ^{t) — ^i(-, t) + i^2(-, t)). By Lemma 17.1 we know that Q is continuous function 
of t > 0, and 

IW)l<^^f^. *>0. (18.2) 
Hence, the Proposition 6.2 is reduced now to the following 

Proposition 18.1. i) Let a function Q{t) e C([0, 00); R^), and *o £ H^. Then the 
equation (18.1) admits a unique solution ^{t) G C([0, 00); if^) with the initial condition 
^(0) = ^0- 

a) If'^Q e Hp and the decay (18.2) holds, the corresponding fields also decay uniformly 
in v: _ 

mm.-,<^^^f$j§i^, t>o, (18.3) 

for \v\ < V with any v e (0, 2m). 

Proof. The statements follow from the Duhamel representation 

*(i) = W{t)^o - [ W{t- s)Q{s) ■ Vp ds, t > 0, (18.4) 
Jo 

where W(t) is the dynamical group (propagator) of the free equation 

i4f{t) = (-A + + iv ■ V)*(t). 
Lemma 18.1. Let \v\ < v with any v e (0, 2m). Then for ^0 e H^. 

WWmoh-p < C{v){l + \t\)-'/^\\^o\U,fs, t > 0. (18.5) 

Proof Note that W{t)'^o = e-^(""^"l^l^/^)*e^^^/2$(t), where $(t) is a solution to free Schro- 
dinger Equation 

i^{t) = -A$(t), $(0) = e*^^/2^o- 

It is wellknown that $(t) satisfies the estimate ||$(i)||i < C(l + |i|)-3/2||$(0)||i,^, t > 
(see for example [15]). □ 

Now (18.3) follows from the condition (18.2), and the Duhamel representation (18.4). 

□ 

19 Appendix 

A. Solitary waves 

Let us to check the last equation of (2.7): 

= J (V^.i(y)pi(y) + V^|J,2{y)p2{y))dy. (19.1) 



Let us transfer to the Fourier representation. Set 

It is easy to compute that 

-ikvipyi + (k^ + rn'^)'ijjy2 ^ -P2, (k^ + m'^)'ijjyi + ikv'ijjy2 ^ -Pi- (19-2) 

Therefore 

J _ -{k^ + m^)pi{k) + ikvp2{k) j _ -ikvpijk) - + m^)p2{k) 
By Parseval identity (19-1) becames 

= y + ^„2P2)^/^ = y (fe2 + ^2)2_(fe,)2 ' 

which is true, since the integrand is odd. 
B. Computing Q{ri,rj) 

Let us to justify the formulas (3.4)-(3.6) for the matrix Q. For j, I = 1,2,3 one has from 
(3.2) and (3.1) 

^i) = diijJv2) - {dji/Jv2, diipvi), (19.4) 

Q(rj+3, T1+3) = {d^jipvi-, dviipv2) - {dv^'ipv2, dy^i^yi), (19.5) 

and 

Q,{Tj, Ti+a) = -{djipyi, dyii)^2) + {dj'ipv2, dyiipyi) + Cj ■ ei. (19.6) 
Differentiating (19.2) we get 

1 _ kjkvTp^i - ikj{k^ + m^)'4)y2 ^ j _ ikjjk'^ + m^)?/'„i + kjkvipy2 • _ -, 9 o 

"^v.y^vl- (^2+^2)2_(fc^)2 ' ^v,^l^v2- (^2 + ^2)2 _ (^^)2 > J " A 

(19.7) 

Then for j, I — 1,2,3 we obtain from (19.4) by the Parseval identity that 

^{Tj,Ti) = j kjki dk{i^yitpy2 - '0«2'0„i) = 0, (19.8) 

since the function ipyc — '4'vi'>Pv2 ~ ''Pv2'>Pvi is odd. Similarly, by (19.5) and (19.7) 
A;,A;,f2z(A;2 + m2)A;^;(|V;,i|2 + |V^„2|')-((A;2 + m2)2 + (A;i;)2)V;^e)(iA; 

(19.9) 



n(r,+3, THS) - - / ((;t2 + ^2)2_(^^)2)2 " 



Finally, by (19.6), 



r kjkiUk^ + m2)(|^^ip + |^„2n +ikvijjyc) dk 

= / ^ (i;. (fa). — + -1 ■ (is'-i") 



Now (3.4) - (3.6) are proved. 



C. Positivity of /i and /2 

Here we check the inequahties which we have used in the proof of Lemma 15.1: 



1 



+ 



1 



M-\v\ki-uj M-\v\ki + uj M-\v\ki 



+ 



( 1- 

\M + \v\ki 



1 



u M+\v\ki+u M+\v\k 



-) >0, 
1^ 



2)A=f i + i 1^) 

\M - \v\ki-uj M-\v\ki + uj M-\v\kJ 

( ^ ^ ^ ^ > fl9 11) 

\M +\v\ki-uj^ M +\v\ki+(jj M+\v\ki)~ ^ ' ^ 



under the conditions \v\ < 2m, < \uj\ < n = w? — First, let us note that the 

expressions in each bracets is positive, since 

112 2a^ 

+ T = 7^ TT, > 



h — a h-\- a h {b + a){b — a)b 

if b—a, b+a > 0, 6 > and it immediately implies that /i > 0. Next, the first summand in 
LHS of (19.11) obviously is not less than the second summand since \v\ki > 0. Therefore 
/2 > and /2 < /i. 
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